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We consider a model in which positive and negative particles with equal densities 
diffuse in an asymmetric, CP invariant way on a ring. The positive particles hop 
clockwise, the negative counter-clockwise and oppositely-charged adjacent particles may 
swap positions. The model depends on two parameters. 

Analytic calculations using quadratic algebras, inhomogeneous solutions of the 
mean-field equations and Monte-Carlo simulations suggest that the model has three 
phases. A pure phase in which one has three pinned blocks of only positive, negative 
particles and vacancies and in which translational invariance is broken. A mixed phase 
in which the current has a linear dependence on one parameter but is independent of 
the other one and of the density of the charged particles. In this phase one has a bump 
and a fluid. The bump (condensate) contains positive and negative particles only, the 
fluid contains charged particles and vacancies uniformly distributed. The mixed phase 
is separated from the disordered phase by a second-order phase-transition which has 
many properties of the Bose-Einstein phase-transition observed in equilibrium. Various 
critical exponents are found. 
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1. Introduction 

While the conceptual framework for the understanding of spontaneous symmetry 
breaking (SSB) occurring in systems at thermal equilibrium is well established, it is much 
poorer for systems far from equilibrium. It was realized recently though, that exclusion 
models with hard core interaction, are simple enough to allow analytic methods, either 
approximative or exact, for the description of their stationary states, yet exhibit a rich 
variety of behaviors even in one dimension. For example SSB of CP symmetry in a 
first-order phase transition was observed in open systems with two [jl], ^], |3| or three 
species 0, |J. 

The aim of the present work is to analyze the structure of the stationary states of a 
one-dimensional three-species exclusion model on a ring with continuous symmetries 
and, in particular, to address the question of the existence of SSB. In this model, 
introduced in |J, positive and negative particles diffuse in an asymmetric, CP invariant 
way. The positive particles hop clockwise, the negative particles counter-clockwise 
and oppositely charged adjacent particles may exchange positions. The model thus 
defined is translationally invariant and the numbers of positive and negative particles 
are conserved. The study presented in || suggests the existence of three phases in the 
stationary state. In particular, as we shall see, for two of these phases charge segregation 
(segregation of species) occurs in the system. Here we confirm these predictions, we 
provide a more complete analysis of these phases and try to unravel the mechanisms by 
which phase transitions occur in such a simple system. A summary of the features of 
the phases and of the phenomenology of the model is given below. We first come back 
to its definition. 

We consider a ring of L sites, labeled by k = 1, . . . , L. Each site is occupied by one 
of the three types of particles, 0, 1 or 2. The particle of type is hereafter referred to 
as a vacancy. The dynamics of the model is defined as follows. At each time step a 
bond is chosen randomly, and the particles a and (3 located at the two ends of the bond 
exchange positions with the rate g aj p, where the indices a, (3 = 0, 1, 2. Therefore, the 
dynamics conserves the number of particles N a of each species. Up to this point, this 
defines a class of models depending on five parameters, since rates are defined up to a 
rescaling of time. Here we restrict the model taking the following rates 
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i.e. the remaining rates go,i an d (72,0 are zero. The model depends on two parameters A 
and q (for convenience, we introduce also the parameter r — 1/q). In order to facilitate 
the intuitive description of the model, we will also denote the particles of type 1 and 2 
by (+) and (— ) respectively and rewrite ( |TTT| ) as 

9+,o = 9o- = A, g- + = 1, g+- = q (1.2) 

In what follows we will consider only the case of equal densities of charged particles that 
we are going to denote by p. With the choice of rates ( |1.2|) , we expect therefore the 
stationary state to be CP invariant and the currents of positive and negative particles 
to be equal. The case of unequal densities is subject of a separate paper @. 

For the sake of clarity, we present here a summary of the main phenomenological 
features of this model, using a simple intuitive presentation. Our knowledge is based on 
the study of the model by exact algebraic methods, mean-field analysis and Monte-Carlo 
simulations, described in the later sections. For a given density p, the phase diagram of 
the model consists of three sectors in the q — A plane. 

The pure phase. (0 < q < 1, any A and p). If one takes a system of size L and start 
from an arbitrary configuration, in the long time regime, the system organizes itself into 
three blocks 

(0)-(0)(+). ..(+)(-)•••(-) (1.3) 

building what we are going to call a bump which covers the whole ring. This may be 
simply understood looking at the rates (|1.2j). Since q < 1, in the average, particles of 
type (— ) will be preferentially located to the right of particles of type (+). Themselves 
will be found to the right of particles of type ( ), while the latter will be found on the 
ring to the right of particles of type (— ). The average waiting time between hops of the 
bumps from one position on the ring to another one increase exponentially with L. Let us 
give an intuitive explanation of this point. Prepare the system as in (|1.3j ). According to 
( |I~2|) the single bond that is allowed to evolve is the one located at the interface between 
the positive and negative blocks. Consider the positive particle at the interface. It will 
begin diffusing through the negative block to its right, but with a bias to the left since 
q < 1, resulting in a restoring force on the particle. As a consequence, the probability 
for this particle to transverse the negative block is exponentially decreasing with the 
size of the system. In the thermodynamic limit, any configuration (|1.3| ) on the ring is an 
absorbing state of the dynamics, i.e. the system does not evolve. Otherwise stated, the 
configuration space is decomposed into a continuous infinity of "ergodic" components, 
i.e. ergodicity is infinitely broken. As a consequence, in the thermodynamic limit, 
translational invariance is broken, the current is zero and the segregation of charges 
maximal, corresponding to a complete charge segregation (hence the name pure phase) . 
Let us observe that, in general, if one denotes by c at p{R, L) any connected two-point 
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function for a system of size L and if instead of the distance R one takes the dimensionless 
variable y = R/L, charge segregation implies that the limits 

lim Ca :l3 (R,L) = Van c a>0 (y,L) = c a> p(y) (1.4) 

L—too L— >oo 

exist (similar limits exist also for any n-point functions). Charge segregation was seen 
in two JI], ||, 0, U and three-species open systems when first-order phase transitions 
take place. In the present case the functions c a ^(y) have a simple expression which 
can be computed using ( |1.3| ). The pure phase was seen also in Ref.[|] for p = 1/3 in a 
different model where the rates are 

01,0 = 00,2 = 02,1 = 1, 00,1 = 02,0 = 01,2 (1.5) 

The finite-size behavior of the present model and that of Ref.|| are different. 

The mixed phase. For fixed A and p, this phase appears as soon as 1 < q and survives 
until q reaches a critical value 

q c = 1 + (1.6) 
y l + 2p K J 

Again segregation of species is observed, though differently from what was observed 
previously for q < 1. In order clarify the physics let us take q larger but close to 1. 
The two blocks of (+) and (— ) particles, which were distinct in the pure phase, merge 
in one region, hereafter called condensate. The rest of the system, hereafter called 
fluid, is mainly occupied by vacancies and traversed by a few (+) and (— ) particles 
distributed in an uniform way. This may be understood as follows. Prepare the system 
in the configuration ( |1.3| ). The only particles that can move are located at the (+)(—) 
interface. The difference is that now the bias, proportional to q — 1, is positive, i.e. 
oriented to the right. Therefore the (+) particle has a finite probability to traverse the 
negative block on its right. As soon as a (+) particle leaves the condensate at its right 
end, it reappears almost instantly at its left end. In other words the positive component 
of the condensate melts on the right but is restored on the left. At the left hand side 
of the condensate the density profile of positive particles is close to one while small 
on the right hand side. The same reasoning holds for the negative particles, we just 
have to exchange the word right with left. Hence it is clear that the condensate as a 
whole moves diffusively on the ring and therefore translational invariance is not broken. 
Finite-size effects are much more important in the mixed phase than in the pure phase. 
In the thermodynamical limit, in the condensate, the densities (+) and (— ) particles 
are both equal to 1/2, hence the condensate can be seen as a two species system with 
open boundaries in the "maximal current phase" . In the "maximal current phase" , the 
current has the expression 

(1-7) 
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which is precisely the value of the current observed in the whole mixed phase. The 
condensate is glued to the fluid where all three species are present and distributed in a 
uniform way. If we let now vary q from 1 to q c , keeping A and p constant, the condensate 
shrinks and the fluid takes over the whole ring at q c . In the whole mixed phase one has 
charge segregation since the limit ( |1.4| ) exists. This brings us to two types of charge 
segregation on the ring. One in which translational invariance is broken (type B), like 
in the pure phase and another one in which translational invariance is unbroken (type 
UB), like in the mixed phase. 

The disordered phase occurs when q > q c , the current has a non trivial dependence 
on the density. The density profiles are uniform and there is no charge segregation. 

A lot of work was done previously on the three species diffusion problem in other 



models, especially in two dimensions (see for example Ref. [1C] and references included). 
The present model is interesting for two reasons. First, this model has the phase 
structure we just mentioned, secondly, the stationary states of this model can be 
studied using the powerful method of quadratic algebras. This not only allows to obtain 
precise results but also allows the introduction of computational methods, like the grand 
canonical ensemble which can not be defined otherwise. This forces us to present our 
results starting with some mathematics. 

In Sec.|| we first review some results obtained in Ref. [[□]] where it was shown that 



for the rates (|1 . 1|) , the probability distribution for the stationary state can be obtained 
by taking traces of monomials of generators of a specific quadratic algebra. We give 
the representations of this algebra. They are in general infinite-dimensional unless A is 
related to q. We next show how to do calculations in the canonical ensemble and explain 
why this approach can be applied only to small lattices. A grand canonical ensemble 
is defined using the algebra. It is stressed that this definition is not unique but very 
useful. We also show a novel application of the algebra which gives the large q (small 
r) expansions of different physical quantities. A connection between the representations 
of the quadratic algebra and the quantum algebra U q so(2, 1) is given in |Appendix A . 



In Sec.|3| we give the values of the currents in the three phases using mainly Monte- 
Carlo simulations for the canonical ensemble and algebraic methods for the grand 
canonical ensemble. We also discuss the problem of fine-size corrections. 

The two-point correlation functions for the pure and mixed phase are presented in 
Sees. | and ||. 

Mean-field calculations are presented in Appendix B| and in Sec.^ first on the lattice 



and then in the continuum. The latter are interesting on their own because we obtain 
stationary solutions of two coupled non-linear differential equations in which the size L 
of the system is a parameter. The mean-field results are compared with various results 
obtained through other methods. 
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In Sec.|7| we define an complex order parameter that we denote by M and compute the 
Ginzburg-Landau potential defined in Ref. || (there it was called free energy functional) 
as a function of \M\. A minimum at a non-zero value of \M\ of the Ginzburg-Landau 
potential would imply charge segregation. This is indeed seen for both the pure and 
mixed phase but not, as expected, in the disordered phase. If we look however for a 
system of size L at the average time T to move from a domain in the complex M- 
plane to another, we observe that T depends exponentially on L in the pure phase 
and algebraically in the mixed phase. This implies that translational invariance is 
spontaneously broken in the pure phase but not in the mixed phase. 

At q c a second-order phase transition takes place. The values of several critical 
exponents are given in Sec.§. It is also shown that, like in the case of Bose-Einstein 
condensation (reviewed in [Appendix Q ), at q c the fugacities do not determine the 
densities. 

A detailed study of the phenomenon of spatial condensation is given in Sec.||. We 
show the L dependence of the compressibility and the Ginzburg-Landau potential as a 
function of the local density of vacancies. We also show that, amazingly, one can define a 
quantity which behaves like the pressure in equilibrium systems. Based on what is done 
in the case of Bose-Einstein condensation, we give a way to eliminate the condensate 
and keep only the fluid in the mixed phase. 

Our conclusions are given in Sec. [Hi 



2. The quadratic algebra and its applications. 

In this section we show how to use the algebraic methods, pioneered by Affleck et al 



12|| for equilibrium problems and by Hakim and Nadal |H| and Derrida et al for 



non-equilibrium problems, to describe the stationary probability distribution obtained 
from the master equation given by the rates ( |1 . 1|) . 

In Ref. it was shown for which rates g a p the unnormalized stationary probability 
distribution P s {/3i, fa, ■ ■ ■ , fa) of a three-state problem with L sites on a ring is given 
by the expression: 

P s (fa, fa,..., fa) = Tr(D 3l Dp 2 ■ ■ ■ D Pl ) (fa = 0, 1, 2) (2.1) 

where the -D Q 's satisfy the quadratic algebra: 

gapDaDp - gpaDpDa = XpD a - X a Df3 (2.2) 



The cases for which the algebra ( |2 . 2| ) exists and has representations with a finite trace 
can be classified depending on how many of the parameters x a are zero. In particular, 
it was shown that for the rates considered in this paper, the probability distribution can 
be obtained by the ansatz Q-Q with x = (See Ref. JTTJ], Eqs.(3.30)-(3.47)). In 
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this case, the quadratic algebra takes the simple form: 

G 1 G 2 =r(G 2 G 1 + X(G l + G 2 )) 
G\Gq = Gq 

GqG 2 = Go (2-3) 

where 

D = d G , D 1 = dtd , D 2 = d 2 G 2 (2.4) 

d Q ,di and d 2 are free parameters (the freedom in choosing x\ and x 2 was taken into 
account). The existence of the free parameters is a consequence of the conservation of 
the numbers of particles N± and N 2 . 



We will first review the representations of the algebra ( [2.3|) obtained in Ref. i I I 
and then explain how they can be used to compute various physical quantities. We will 
also show how to use the algebra ( |2.3j ) to obtain a large q (small r) expansion and derive 
expressions for the current. It turns out that in some cases, G\ and G 2 can be expressed 
in terms of generators of the quantum algebra U g (so(2, 1)). In Appendix A we show 
how this information can be used. 



2.1. Representations of the quadratic algebra. 

As discussed already in Ref.[|ll|], we are interested for a given value of A and r in 
representations of the smallest dimension. 

A one-dimensional representation (this is equivalent to a product measure 
probability distribution) is obtained for any A and r: 

G = 0, G x = 7i, G 2 = l2 (2.5) 

with 

7i(l — r) — Ar 

Here 71 and 72 are c-numbers. Since Go — 0, one has no vacancies and one has only an 
asymmetric exclusion process. 

If Go is non-vanishing it is convenient to denote 

G 1 = A + U, G 2 = A + V (2.7) 

and a representation is given by: 



Gn 



/ 1 







V ; 



A 



/ si 
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a 4 
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( t x • • • \ 
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t 3 
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si 

s 2 

s 3 



:2.81 



with 

a k = r[(X + l){k-l} r -{k-2} r ] 

s k t k = {k} r [(2A + l)r - 1 + (r(A + 1) - l) 2 {A; - l} r 

We have used the notation 

{k}r ' 



r — 1 



(2.9) 



(2-10) 



The representations given by Eq. (|2.8| ) are not necessarily infinite dimensional. One 
can obtain an n-dimensional representation by finding the solutions of the equation: 



Sn tn . 

Using Eqs. (pO|) and ( |2.11| ) we obtain: 

1 — r 



A 



(2.11) 



(2-12) 



r(l + r 2 ) 

The n-dimensional representation has the following matrix elements (k = 1, 2, . . . , n) 
1 + r 2 



1 

1 + r 2 

(1 — r fc ) (1 — r k ~ n ) 



1-n \ 2 
7" 2 



(2-13) 



Let us observe that from Eq.( |2.12| ) we get a one-dimensional representation, n — 1, with 
Gq = "1 = ^2 = 1 if 



q = 2A + 1 



(2.14) 



Notice that in this is case Go is non-zero! A two-dimensional representation, n = 2, is 
obtained if 



g=(l + A) 2 
and one can choose the basis such that one has 



A 



1 

xff 



, u 



1-Vr 




V 




r - 1 



(2.15) 



(2.16) 
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A second two-dimensional representation is also obtained for any A if r = 0. In this case 
one has to take r = in the expressions for A, U and V in Eq.( [2.16 ). This observation 
will be useful in SeaEO. 



To sum up, we have a one-dimensional representation with Go = for any A and 
q and a second one-dimensional representation with Go 7^ when q and A are related 
by Eq. Q2.14l ). We also have two two-dimensional representations. One corresponds to 



Eq.( |2.15D , the second one is obtained for r = 0. For n > 3, one has only one n- 
dimensional representation when Eq.(|2.12D is satisfied. For generic q and A and Gq 7^ 



the representations are infinite-dimensional. 

Due to the special structure of the representations with trace of the algebra (|2.3|) , 
the expression ( ]2.1|) has to be used with care. If one looks at configurations with no 
vacancies, one has to use the one-dimensional representation (|2.5| ). For any configuration 
which contains at least one vacancy one has to use the representation given by Eqs. (|2.7|) - 
( p.8|) which, as discussed, can be finite or infinite-dimensional depending on the values 
of q and A. 

2.2. The canonical ensemble. 

In this case the numbers of particles N\ and N2 are given. The normalized probability 
distribution can be obtained from Eqs.(|2.1|) and ( |2.4| ): 



-Pvi,V 2 (A'^2, • • • ,Pl) - Nl >JV2 I i N 1 t N 2 (Pl,02, ■ ■ ■ , (2-17) 

where Z^ 1 ' 2 is a normalization factor and 



-^Vi,7V 2 (^l> $2, • • • , Pl) — ^{Gpi Gf3 2 ■ ■ ■ Gp L ) 

x5(E 5(P k - 1) - N{) 5C£ 5{(3 k - 2) - N 2 ) (2.18) 

k=l k=l 

where 5{x) is the Kronecker function. 

The currents are going to play an important role in this paper. They are defined 
as follows: The current for the particles of type 1, taken for convenience on the link 
between the sites L — 1 and L, is: 

h = A 5(J3 L -! - 1) 5(J3 L ) +q5(p L - 1 - 1) 5(J3 L - 2) - 5([3 L - 1 - 2) 5(J3 L - 1) (2.19) 
The current for the type 2 particles has a similar expression. Using the Eqs. (|2.17| )- (|2.19| ) 



and (|2.3j ), the expression for the average current is: 



x5(J2S(Pk ~ 1) - #1 + 1) 6(£ 6{/3 k - 2 ) - ^2) 

k=l k=l 
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+Tr(G / g 1 Gp 2 ■ ■ ■ Gp L _ 2 (Gi + G 2 )) 

x8(£6{(3 k - 1) - Nt + 1) 6(E 6(/3 k - 2) - iV 2 + 1)] (2.20) 
fc=i fc=i 

It is useful to employ the Fourier transform of the Kronecker delta function 

1 L_1 2tt\ 

6(a) = - 52 exp(-j-ka) (a = 0, . . . , L - 1) (2.21) 
L k=o L 

and to define 

B% =G + exp(^r)G 1 + exp(^ S )G 2 . (2.22) 

In this way one can get a compact expression for the current of type 1 particles. The 
expressions for the type 2 particles and for the normalisation factor Z^ 1,N2 can be 
obtained in a similar way. One obtains: 

= i J>xp (-27ri (r§ + S §)) Tr ((B$f) ( 2 . 23 ) 
* = ^ |>P (-2.1 (r^l + ^)) Tr ((^-2))-^)) (2.24) 

^ = ^ |>P (-2- (rJ^ + s^l)) Tr ((Sg^)^^) ( 2 ' 25 ) 

Analogous expressions can be written for the correlation functions. 

In a concrete calculation (q and A are given) we take G , G\ and G 2 from Eqs.fl2.8[) 
and fl2.9| ). This gives the matrices B^ L ) of Eq.( |2.22| ) and we are left with the calculation 
of the traces appearing in Eqs. (|2.23|) - (|2.25|) and the double summations over r and 



s. We have been unable to do any analytical calculation in this way in the general 
case (for finite dimensional representations this is probably possible). This takes us to 
numerical calculations. If we take Ni + N 2 ^ L, the trace operation in Eq.( 2.20| ) is 



well defined and any divergence in Eqs. (|2.23| )- (|2.25| ) has to cancel, this does not make 
however numerical evaluations easy. To circumvent the problem of computing traces of 
infinite dimensional matrices we will assume that we have at least one vacancy and take 
in Eq.(^23|) Tr^o^' 1 ^ -1 ) = {(B^ 1 ))^ 1 }^ instead of Tt((B^J) l ) and L - 1 
instead of L. The expressions for the currents have to changed in the same sense. Due 
to the fact that G is a projector (see Eq. ( [2.8|) ) this procedure has also the advantage 
that for a problem with L sites, only the matrix elements with k < L (see Eq. ( j2.8| )) 
appear in the evaluation of the traces. In this way we have been able to compute the 
current up to L = 30 with high precision. This allows in principle, using extrapolants 
I4| , to obtain the current in the thermodynamical limit. As will be shown in Sec^Jthe 



current might converge very slowly towards its thermodynamical limit and therefore in 



11 



these cases this method is useless. As a matter of fact, throughout this paper, we have 
used Monte-Carlo simulations for the canonical ensemble. The precision is much worse 
than using the analytical method but one can reach much larger values of L. 



2.3. The grand canonical ensemble. 

As opposed to equilibrium states, the definition of a grand canonical ensemble for 
stationary states is not unique. In principle, for the master equation one can give, 
at t = 0, an initial probability distribution in which the number of particles N% and N 2 
are not fixed but the "chemical potentials" \i\ and fi 2 are. Such an initial probability 
distribution has the expression: 

P%,jPi,---,PL,t = 0)= £ QWi,...,(3 L ;N 1 ,N 2 )e^+^ (2.26) 

N lt N 2 

where Q is an arbitrary function. This initial condition determines a certain stationary 
probability that can be considered as a grand canonical ensemble. The N\ and N 2 
dependence of the functions Q remains however arbitrary even if we make sure that 
the average densities p\ =< Ni > /L and p 2 =< N 2 > /L are properly obtained. We 
will adopt a different approach looking directly at the canonical stationary distribution 
given by Eq. ( |2.17| ) and consider the grand canonical ensemble defined by the expression: 

.".fc) = ir E ^y'-'^ ^^ (2.27) 

N lt N 2 JNi,N 2 

where we have used Eq.( j2.1§ ). In Eq. ( 2.271) the factors /jvi V2 are arbitrary. We will 
take from now on: 

A similar definition was used in jl5fl . This gives: 



P 9 ^ 2 (Pi, ■■■,Pl) = ^Tr(^ • • • Ep L ) (2.29) 



where 



and 



E = G , E 1 = z 1 G 1 , E 2 = z 2 G 2 (2.30) 



Zl = e^, z 2 = e" 2 (2.31) 

are the fugacities. As opposed to the canonical ensemble, the grand canonical partition 
function has a simple expression: 

z ™ = Tr ^ ^ 32 ^ 
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where 

C = E + E 1 + E 2 . (2.33) 

That this choice for the definition of the grand canonical ensemble is not unique can 
be seen in a different way. Instead of using the D's of Eq.( |2.4| ) with d = di = d 2 = 1 
(for the canonical ensemble this choice is not important), we could have taken d = 1 
(one of the di can be chosen at will) but keep d\ and d 2 as free parameters. In this case, 
Eq.( |2.30D is replaced by 



E = G , E 1 = d 1 z 1 G 1 , E 2 = d 2 z 2 G 2 (2.34) 

and all the results for finite-size systems will be affected, but hopefully not the results 
for the thermodynamical limit. We want to stress this point because the interpretation 
of finite-size corrections for the grand canonical ensemble has still to be understood. We 
have arbitrarily taken the choice given by Eq. (|2.28|) . Our choice promotes the algebra 
( p. 3D to a special role, since the isomorphic algebras obtained using various values of 
the di in Eq.( [2.4|) give different definitions of the grand canonical ensemble. We have no 
good reason to prefer one algebra to the other except the fact that the representation 
of the algebra (|2.3|) can be connected "naturally" to known representations of quantum 
algebras (see Appendix A). 

The expression ( |2.29|) reminds us of an analogous expression in equilibrium statistical 
physics and along this line of thought, one could perhaps find a better justification for a 
proper definition of the grand canonical ensemble. That this idea is not too far fetched 
will be seen in Sec.|9]. 

We now proceed with the description of the grand canonical ensemble. As usual the 
fugacities are fixed by the densities: 

» = — = z^kf- (2 - 35) 

The fluctuations of the densities are: 

< (AiV,) 2 >= dl °* Z \ (2.36) 
dlogZi 

Using Eqs.( |2.29| ) and (|2.33|) , one gets the two-point functions: 



CaAR) = -^TriGaC^G^- 11 - 1 ) (2.37) 

Notice that if the matrix C has the two largest eigenvalues Ai > A2 and if these 

eigenvalues are not degenerate, at large distances R the two point function has the 
expression: 

c a>l3 {R) oc exp(-Mi?) (2.38) 
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where the inverse correlation length M is: 

M = log^ (2.39) 



This implies that for finite-dimensional representations (see Eq.( 2.13f) ) one expects finite 



correlation lengths. This conjecture can be proven rigorously but it is a lengthy exercise 
that we have not done. The expression of the currents can be derived using Eqs. (|2.19Q 



and (2J). We get 



Ji = ^SuTr ((Go + z 2 (G 1 + G 2 ))C L - 2 ) . (2.40) 



The expression for J 2 is obtained by permuting 1 with 2 in Eq.( j2.40| ). If we consider 
the case in which the densities of the particles 1 and 2 are the same, the expression of 
the currents takes a very simple form: 

Ji = J2 = A^fiJ (2.41) 



where 



e M . (2.42) 



This makes the calculation of the current much simpler in the grand canonical ensemble. 

In the expressions written above (see ( |2.32p , (|2.37| ) and ( |2.4(J| )), we have made the 
assumption that the trace operation exists. This is clearly the case if one has finite- 
dimensional representations. The case of infinite-dimensional representations has to be 
considered with care. In the grand canonical ensemble, for a given number of sites L one 
sums over the values Ni and N 2 . In this sum there is also the term with N% + N% — L 
where no Go matrix appears (|2.18|) and one has only products of G\ and G 2 matrices 



given by Eq. fl2.8Q which have no trace. One has to use for this case the representation 
(2.5) which makes the calculation cumbersome. A simpler way is to define the grand 



canonical ensemble such that one has at least one vacancy. That is what we have done. 
In this case, one has to sum over all possibilities for the relative position of the Go matrix 
and the operator whose expectation value one wants to compute (see Eqs. (|2.32 ),( 2.37|) 



and ( |2.40| )). For each term, taking into account the form of the matrix Go (see Eq. fl2.S| )), 
one replaces the trace operation by taking an appropriate (1,1) matrix element. Like in 
the canonical ensemble, we have not been able to do analytical calculations and had to 
use the computer but in the grand canonical ensemble we were able to reach L = 300, 
as opposed to L = 30 for the canonical ensemble. 

2.4- The large q expansion in the grand canonical ensemble. 

A useful application of the algebra ( |2.3| ) is that one can do a large q (small r) expansion 
in an algebraic way. We will show how to obtain the physical quantities up to order r. 
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We first note that, if we keep terms up to order r in the matrix elements ( |2.9| ) which 
give Gi and G 2 , the only non- vanishing elements are 



ai = 

Siti 



(2A 



a 2 = (l 
IV - 1 



A)r 



dfc = Xr (k > 3) 



S2h 



(2.43) 



Since one is interested in the smallest representation, it is enough to take a 3 x 3 
dimensional one. One can check that taking this representation one obtains: 



GiG 2 - r (G 2 Gi + X(Gi + G 2 )) 



f 2 



(0 








(2.44) 





-(1 + A) 
V (1 + A) 

which shows the consistency of the approach since the right-hand-side of Eq. (|2.44j ) is 
of order r 3//2 . We have chosen a representation in which U = —V T (see Eq.( |2.8| )). Notice 
that if we neglect the terms of order r, we obtain the two dimensional representation 
( j2.16| ) in which we take r = 0. 

We will consider only the grand canonical ensemble and assume that the densities 
of particles 1 and 2 are the same: 



P = Pi = P2 

The matrix C given by Eq.( |2.32 ) becomes 



C 



( 



\ 



l + Zl + z 2 Zl {l-^r) 
-z 2 {l-^ 1 r) (1 + A)(*i + z 2 )r z x A 

— Zir? X(zi + z 2 )r J 



\ 



(2.45) 



(2.46) 



We are interested to obtain the coefficients of the r-expansion defined by the 
expressions: 

Xi(z u z 2 ) = Xf\z u z 2 ) +rX { i\z l7 z 2 ) 



Zl/2 



M 

'1/2 



(P) 



rz 



(i) 
1/2 



(P) 



J l = jf\p)+rJ\'\p) 



(2.47) 



where Xi(zi,z 2 ) is the largest eigenvalue of C. We didn't look for higher order terms, 
which would imply taking larger matrices than 3x3. 

We start with the leading order in which case we have to deal with a 2 x 2 matrix. 
The two eigenvalues of the matrix C are: 



\(0) 
A l/2 



2 ( 1 : J : 2± 



'1 + Z! + z 2 f - Az lZ2 



Using Eqs. (l2T32D and flEggp we get: 



'1/2 



(P) 



;i - 2 P ) s 



(2.48) 



(2.49) 
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Introducing the expression ( |2.49| ) in Eq.( |2.48D we find the eigenvalues of the matrix C 
at a given density: 



\f{z,z) 



1 - 2p / 



(2.50) 



Using now Eq.( 2.39|) , we derive in zero-th order the expression of the inverse correlation 
length: 

/ l-p\ 



M = 2 log 



P 



(2.51) 



It is interesting to observe that the "mass" vanishes for the maximum value of the 
density (p = 1/2) and diverges at small densities. We can now use the expression of the 
partition function 



(0) 



and Eq. (|2.41 ) to get the expression of the current: 



Ji 



(0) 



_Xp_ 

1-p 



(2.52) 



(2.53) 



To get the next order terms in r, we use the expression (|2.46 ) for the matrix C and get 
the correction to the largest eigenvalue. We obtain: 

(1 , ^_ a(AS°)) 2 -/?AS°) + 7 



3(Ai 0) ) -25\f ) + z 1 z 2 



(2.54) 



where 



a = (1 + 2A)(5- 1) 

P = (l + 2X)5(5 - 1) - 2\ Zl z 2 

7 = ((1 + X)5 - \)z lZ2 

5 = 1 + z 1 + z 2 

It is obvious how to proceed further. We get: 
(1) _ l + 2A-2p-2Ap 

Al {z ' z) ~ p (T^2pJ 

(1) _ A + (l-3p + p 2 )(l + A) 
z — p 



(2.55) 



1 - pf (2p2 + 1 - 2p) 



(2.56) 



(i> 



2 1 + 2A - p(6 + 9A) + p 2 (12 + 13A) - p 3 (6 - 2A) - p 4 (6 + 8A) + 4p 5 (l + A) 
" P (l-p) 4 (2p 2 + l-2p) 

One interesting result that one can read from Eq. ( |2.56| ) is that, at fixed A, the current 
reaches the asymptotic value for q —* oo from below. 
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Figure 1. The current J as a function of q for A = 1 and two densities (p = 0.2 (+) 
and p = 0.4 (o)) The mean- field values are given by the dashed lines (a) (p = 0.2) 
and (b) (p — 0.4). The line (c) is given by Eq^JiJj). The data are the results of 
extrapolations using the grand canonical ensemble and L up to 200. 



3. The currents. 

From now on, the particles of type 1 and 2 will be called positive respectively negative 
because we think that it is easier to visualise the physics in this way. Their densities 
will be denoted by p respectively m (the density of vacancies will be v ). If the densities 
of positive and negative particles are equal (this is going to be mostly the case), we will 
write 

p = m = p . (3.1) 

It is useful to give the expression of the currents in the mean-field approximation. 
From Eq.( p.l9| ) we get: 

\pv + (q — l)pm 



J, 



Xmv + (q — l)pm 



where obviously 



p + m + v = 1 



(3.2) 



(3.3) 



As is shown in [Appendix B| (see also Sec.^]), the mean- field equations have not only 
homogeneous solutions. Here we refer to the homogeneous solutions. We will give first 
some pictures showing the q dependence of the currents at fixed values of A and densities. 
These pictures are just a few from the data we collected. Our conclusions will be based 
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on all the available information. The thermodynamical limit for the values of the current 
have been obtained using the algebraic approach to the grand canonical ensemble (see 
Sec. |2.3| ) up to L = 200. The large L limit has been obtained using extrapolants ||14|| . 
All the results have been checked using Monte-Carlo simulations (this corresponds to 
the canonical ensemble). 

Since we consider the case in which the densities of positive and negative particles 
are the same and therefore one has: 

J + = J_ = J (3.4) 

In Fig.[l| we show the current for fixed A = 1 and two values of the density (p = 0.2 and 
0.4). The large q behavior of the current is not shown in the picture since its expression 
is known from Eq. ( [2.471 ). We distinguish three domains that for good reasons we are 
going to call phases. The names for the phases will look natural in the next sections. 

a) Pure phase (0 < q < 1) where the current is zero. 

b) Mixed phase (1 < q < q c ) where the current is given by 

J="-^ (3.5) 

(see the line (c) in Fig.|l]), for both densities. The critical value of q denoted by q c 
is density dependent. 

c) Disordered phase (q > q c ) where the current is density dependent. For large values 
of q, the current is given by the large (/-expansion (see Eqs.( j2.47| ), ( |2.53| ) and ( |2.56| )). 



The solid lines (a) and (b) in Fig.^Jare the mean- field values obtained from Eq.( |3.2| ). One 
can notice two facts. First, for q = 3 mean-field gives the current correctly. This should 
not come as a surprise since from Eq.(|2.14|) with A = 1, we learn that the quadratic 



algebra has a one-dimensional representation, which implies that mean- field is exact. 
Next, we observe that, for both densities, mean- field gives correctly the currents for 
values of q close to q c . This result has deeper consequences as will be seen in Sec.||, where 
it will be shown that at q c the current is given by mean-field and that the discrepancy 
shown in Fig.l comes probably from finite-size effects (one extrapolates from L = 200). 

Having in mind Eq.( [3.5|) , in Fig.^| we show the quantity 4 J / (q — 1) as a function q 
for fixed density (p = 0.2) and three values of A. The same three phases are present. 
We notice that the current vanishes in the pure phase, that it is given by Eq.( |3.5| ) in 
the mixed phase and that q c increases with A. For A = the current vanishes for any 
value of q (see Eq.(|2~20|)). 

Similar calculations for other densities and values of A have shown that in the 
pure phase the current vanishes and that in the mixed phase it is given by the simple 
expression ( |3.5| ) independent of density and A. This is a very unexpected result. The 
values of q c are dependent on A and on the density. In Sec|] we will derive the expressions 
of q c (see Eq. fl6.26p for equal densities) and the finite-size corrections. The interested 
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Figure 2. The q dependence of 4J/(g — 1) for fixed p = 0.2 and three values of A 
(A = 0.5, 1 and 2). The data are obtained like in Fig. 1. 
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Figure 3. Large L limit of the connected nearest neighbour correlations c+,o(l) = 
Co,-(l), c ,+ (l) = c_ j0 (l), c_ + (l), c + ,+ (l) = c__(l), c+_(l) and c ,o(l) for A = 1 
and p = 0.2. The data are the results of extrapolations using the grand canonical 
ensemble and L up to 150. 



reader can already have a look at Fig.0 where the density dependence of the current is 
shown for various lattice sizes. 

It is obvious from the Figs.|I| and ^| that at q c the derivative of the current has 
a discontinuity. In Fig.^| we show other local quantities, which are various connected 
correlations functions at distance 1 (one lattice site) for A = 1 and p = 0.2. The 
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large L limit was obtained using the grand canonical ensemble extrapolating from data 
obtained up to L = 150. One notices that the derivatives of all these quantities have also 
a discontinuity at q c = 1.62 (see Table [I]). The Fig.|3| shows also an interesting fact: for 
reasons which remain a mystery, the connected function c_ + (l) practically vanishes in 
the disordered phase such that the discontinuity of its derivative is easier to determine. 
This allows the best estimate of q c = 1.62 ± 0.04. 

In the next sections we are going to discuss at length the three phases. It is however 
very instructive to describe what one can see on the screen of the terminal while looking 
at the succession of configurations of a single Monte-Carlo run. This is of course a 
qualitative picture only. First, let us consider the pure phase. In order to fix the 
ideas, assume q = 0.8, A = 1 and p = 0.2 (see Fig.[l]). For lattices lengths L < L\ 
(Li « 20), one sees a uniform distribution of positive, negative particles and vacancies. 
For L\ < L < L2 (L2 ~ 30) one observes a separation into two domains: one containing 
only positive and negative particles and one in which one has predominantly vacancies 
but also some charged particles. For L > L 2 one sees three domains one purely positive, 
one purely negative and one purely neutral (hence the name pure phase) . These domains 
look like pinned in space (they move very slowly from one place on the ring to another 
one). In the mixed phase, (take q = 1.2, A = 1, p — 0.2, see again Fig.[l]) the evolution 
with L of the relevant configurations is similar but by no means identical. Again for 
L < L\ [L\ m 80) the distributions are uniform, For L\ < L there is again a separation 
into two domains. In the first, one has only positive and negative particles with non- 
uniform distributions, we call it condensate. In the second domain, positive, negative 
particles and vacancies are uniformly distributed, we call this domain fluid. Increasing 
L, the concentrations inside the condensate change slowly and get more uniform and 
the fluid survives as a domain with a uniform distribution of the three species. This is 
one major difference between the pure and mixed phases. A second one comes from the 
observation that these profiles are not anymore pinned in the space but have a Brownian 
motion on the ring. 

Before closing this section we would like to discuss the problem of finite-size 
corrections to the current. 

It is amusing to start with a simple problem. Let us take A = 1 and q = 3 and 
equal densities p. As already observed before, in this case one has a one-dimensional 
representation of the quadratic algebra ( |2.3| ) which can be obtained using Eq. ( |2.13| ): 

G = Gi = G 2 = 1 (3.6) 

One can then do the calculations in the canonical ensemble where one gets in leading 
order 

J = p(l + i) (3.7) 
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Figure 4. The deviation of the current J(L) for L sites from its asymptotic value 0.05 
in the canonical ensemble (q — 1.2, A = 1, p = 0.2). The values obtained using Monte- 
Carlo simulations (+) are shown together with those using the algebraic approach (o, 
up to L = 30) and with those from inhomogeneous mean-field solutions (dotted line) 



as described in Appendix B 



and in the grand canonical ensemble where one obtains 
J = p 



(3. 



for any number of sites. This result is surprising since the parameters are such that we 
are in the disordered phase and one would expect exponential corrections (as will be 
shown in Sec.|8| the correlation length is finite in this phase). 

We have looked in several cases to the finite-size correction terms and the results 
are given below. All data are for A = 1, p — 0.2. First we consider the grand canonical 
ensemble: 



q 

Q 
Q 
Q 



0.8 
1.2 
2.2 
3.8 



J = 0.09exp(-0.055L) 
J = 0.05 + 1.6/L 
J = 0.18158 + 0.04/L 2 
J = 0.21065 — 0.016/L 2 



(pure phase) 
(mixed phase) 
(disordered phase) 
(disordered phase) 



(3.9) 



For the canonical ensemble, it is harder to estimate the correction terms. In order 
to illustrate the problem, we take q = 1.2 in the mixed phase. In Fig.^ we show the 
deviation of the current from its asymptotic value (0.05) as a function of the number of 
sites L. The Monte-Carlo data are shown together with the results obtained using the 
approach presented in Sec. |2.2j . Looking at the Monte-Carlo data first (up to L = 200) 
one observes that one gets large corrections. (The Monte-Carlo data are not precise 
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Figure 5. The y dependence of the two-point correlation functions co,o (a), c + ._ 
(6) in the pure phase (q = 0.5, A = 1, p = 0.2) and lattice sizes L = 75(o, □) and 
L = 100(+, x). (Monte-Carlo simulations). 



enough to use extrapolants). On the other hand, the data obtained in the algebraic 
approach are in a range of L in which any extrapolation is bound to give wrong results. In 
Fig}4| we also show the values of the current obtained from the inhomogeneous solutions 
of the mean-field equations (see Appendix B). At large values of L the current obtained 
in this way is in agreement with the Monte-Carlo data, but at small values of L, they 
show a curious behavior. The explanation of this phenomenon and a more extensive 
discussion of the finite-size corrections in the mixed phase for the canonical ensemble 
can be found in Sec.^|. There is a case in which the finite-size corrections converge 
very nicely. This is the case q = 1. Moreover, in Sec.|6] we will even show how to get 
them theoretically. Here, we just give the Monte-Carlo results, again for A = 1 but for 
different densities: 



p 


= 0.1 


J = 


5/L 


p 


= 0.2 


J = 


2.5/ L 


p 


= 0.4 


J = 


1.25/ L 



The reader can easily guess the analytical expression of the correction term. In order 
to help the reader, we have not given the errors in Eq. ( |3.1C| ) (they are small). 
We will now discuss in detail the three phases. 
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4. The pure phase (0 < q < 1) 

In the last section we have mentioned that in the pure phase the current vanishes 
exponentially with the lattice size for any values of the densities (they might be different) 
and of A. We will now get a complete picture of this phase (it is very simple). 

At q = 0, a single vacancy is sufficient to break the translational invariance of the 
system. At a finite density of vacancies, the ground state of the time evolution hamilton 
operator is infinitely degenerate for L = oo, each configuration of the kind 

(0). (0)(+) •••(+)(-) ... (-)(())■ ■■(())(+) •••(+)(-) •••(-)••• (4.1) 

being a stationary state. At q different from zero and finite L, if one starts with an 
arbitrary configuration, the system organizes itself into only three blocks 

(0)-(0)(+). ••(+)(-)•••(-) (4.2) 

which cover the lattice. Translational invariance is respected but as shown in Sec.^| the 
blocks hop from one position to another with a flip time which increases exponentially 
with L. This implies that for q < 1, translational invariance is spontaneously broken. 
Since each of the three blocks contains one kind of particle only, we call it the pure 
phase. 

We postpone the discussion of the flip times and of the behavior of the order- 
parameter which describes the breaking of translational invariance. This is going to be 
done in Sec.|7] in conjunction with a similar problem but with different results for the 
mixed phase. 

In Fig.|5| we show the two-point correlation functions co,o and c +i _ defined by 
Eq. ( [2.371) n °t as a function of the distance R measured in lattice sites but of the 
macroscopic distance 

y = R/L (4.3) 

For any fixed y, the data converge exponentially to their thermodynamic limit. This 
explains why the correlations functions for so small lattices [L = 75 and 100) coincide 
already. The expression of the correlation functions for very small lattices can be derived 
from the mean-field calculations presented in Sec.§. It was checked using Monte-Carlo 
simulations that these calculations give the correct picture of finite-size effects. 

We invite the reader to compute the two-point functions using the three blocks 
picture ( f4.2|) and he will discover that he will get precisely the results presented in 
Fig.|5]. This applies also to the other two-point functions. Obviously, the correlation 
functions dependent only on the densities but not on A. 

That in the pure phase one has charge segregation is obvious. Not only the 
correlation functions do depend on y in agreement with the definition given in SecfJ, 
but the configurations described by ( |4.2| ) give the simplest example of segregation. The 
mixed phase will provide another example. 
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5. The mixed phase (1 < q < q c ) 
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Figure 6. The correlation functions co,o (a) and c +! _ (b) obtained from Monte-Carlo 
simulations for p = 0.2, q = 1.2, A = 1, and L = 400. The solid lines are obtained 
using the mean-field results presented in Sec.^. 

In Sec.|3| we have seen that in the mixed phase, for any p and A, the current is given 
by the expression ( |3.5| ). We are going to analyse this phase in more detail. In Figure ^| 
we show for q = 1.2, A = 1 and p = 0.2 (one can check from Fig.0 that we are in the 
mixed phase) the correlation functions Co 5 o and function of y for L = 400. The 

"macroscopic" variable y is defined in Eq. (|4.3| ). The correlation functions were obtained 
using Monte-Carlo simulations. As will be explained in Sec.^] one can't use the grand 
canonical ensemble to compute correlation functions for large distances (this is not the 
case for the current) in the mixed or pure phases but only in the disordered phase. 

It is interesting to compare the correlation functions in the mixed phase and in the 
pure phase (see Fig.|5|). They look similar. The co,o function has the same shape in both 
cases although the values of co,o are different. For example the widths of the plateau 
have different values. The c +i _ functions have a slightly different shape (compare the 
behaviors for y < 0.4 and y > 0.7). One could be tempted to conclude that there are no 
major differences between the two phases. This is not the case. As we mentioned, in the 
pure phase the data presented in Fig.|5] represent already the thermodynamical limit. 
The data presented in Fig.|6| although obtained for a large lattice do not represent the 
thermodynamical limit. There is a slow change with L of the behavior of the correlation 
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functions not shown in the figures. We took data for lattice sizes up to L = 800, the 
L dependence is well taken into account by the mean-field calculations presented in the 
next section. That these mean-field calculations make sense can already be seen in Fig.|6] 
where the solid curves which are on top of the Monte-Carlo data are mean-field results. 
After presenting the mean-field results we are going to compare them with other data 
obtained from Monte-Carlo simulations (see Figs.^] and |9]). 

At this point the reader should accept the fact that the correlation functions converge 
at large L to non-trivial functions of y (this is going to be explained in the next section) 
and therefore in the mixed phase, like in the pure phase we have charge segregation. 

As we are going to prove in Sec.[7|the major difference between the two phases is that, 
as opposed to the pure phase, in the mixed phase translational invariance symmetry is 
not spontaneously broken. 



6. Inhomogeneous solutions of the mean-field equations. The pure and 
mixed phases. 

We have postponed the presentation of mean-field for the following reason. We had first 
to present the empirical data for the correlation functions in order to know how the 
mean-field calculations have to be done. Since for a given lattice size L, the correlation 
functions depend smoothly on the "macroscopic" variable y (0 < y < 1), this suggests 
to work in the continuum taking y as variable and L as a parameter in the mean-field 
equations. The previous mean-field calculations in similar problems (see for example 
Refs. [[TIJ and JT7J) were done in a different way. The mean- field equations for the 



stationary state are a couple of non-linear differential equations and it is not clear a 
priori how to solve them. The inspiration came from getting, numerically, the solutions 
on small lattices. The results of these investigations as well as the mean-field equations 



in the continuum are presented in Appendix B| . About the relevance of mean-field to 



our problem, we are going to comment later. 

In [Appendix B| it is shown that the for given values of q, A and densities in the pure 



or mixed phases the mean-field equations on the lattice have besides an homogeneous 
solution (which was used in Eq. ( |3.2|) ) inhomogeneous ones. These solutions are "pinned" 
in space even for a finite lattice. Otherwise speaking, in mean-field translational 
invariance is broken in both the pure and mixed phases. These conclusions are based on 
numerical results. One remarkable property of the inhomogeneous solution (see Figs. BT 



and |B^) is that it shows two domains one with charged particles only (we called it 
condensate) and one with charged particles and vacancies (we called it fluid). In this 
section we are going to present the stationary solutions of the mean-field equations ( B.3 



B.6| ) which we were able to obtain by making the Ansatz of the existence of a condensate 



and a fluid. A stationary non-homogeneous solution of the mean-field equations will be 
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contain at least one bump. In the present paper we consider the solutions with one bump 
only since they are relevant for our problem. Solutions with more bumps are shortly 
mentioned in |Appendix B| and discussed in detail in [|7| . These solutions always breaks 



translational invariance. We will assume that mean-field gives the correct description 
of our physical problem (we will return to this point at the end of this section) if its 
solutions are used in different ways for the pure and mixed phases. 

We will work with "macroscopic" variables like in equation (O). In this case the 



ring has a perimeter of length 1. For the mixed phase we will assume that with equal 
probability for each point z on the ring (0 < z < 1) one has a bump which contains a 
condensate of length a and a fluid of length b (this implies that one has translational 
invariance) and that mean-field describes both of them. For the pure phase we take 
the results of mean- field at face value: the bump is "pinned" in space and translational 
invariance is broken. In both phases the values of a and b are L-dependent and their 
values in the thermodynamical limit will be denoted by A and B. Obviously 

a + b = l. (6.1) 

It is convenient to choose a coordinate y with the origin at the beginning of the 
condensate. Thus, if < y < a we are in the condensate, and if a < y < 1 we 
are in the fluid. The names condensate and fluid are naturally related to Bose-Einstein 
condensation. As will be shown in Sec.|9] a similar phenomenon takes place in our system. 
We assume that there are no correlations in the condensate and the fluid so that we can 
apply mean-field theory. The present calculation will give not only the concentrations of 
particles but also the values of the current and the length b of the fluid as a function of L. 
We start with the fluid and assume that on the segment of length b the concentrations 
of positive and negative particles, p+(L) and p^(£), and implicitly of vacancies, Pq(L), 
are independent of y but L dependent. One has: 

p«(L) + p fl (L) + ^(L) = l (6.2) 

In the condensate, one has no vacancies, the concentrations of positive and negative 
particles (p+(?/,L) and p c °(y,L)) are y and L dependent and satisfy the relation: 

p c °(y,L)+p c _°(y,L) = l (6.3) 

We will consider three cases. We take the same density of positive and negative 
particles 

1 - v 

P = (6-4) 

and q greater than 1, next q = 1 and q < 1. Various predictions obtained from mean-field 
will be compared with the results obtained from Monte-Carlo simulations. 
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Figure 7. The current as a function of the density. The data are obtained from 
Monte-Carlo simulations for q = 1.2 and A = 1 and L = 100 (top), L — 200 and 
L = 400 (bottom). The horizontal line gives the value of the current in the mixed 



phase (see Eq. (|3.5| )) 



6.1. q>l (mixed phase) 

6.1.1. Profiles for finite L In the fluid we have 



pl = p a = \a- P a o) 



and 



bpl = v 

The current is 

J = Ap^ + (g-l)p> fl 
For later convenience, it is useful to define 



a 



/ 4J 
q-1 



(6.5) 



(6.6) 



(6.7) 



(6. 



It is clear that if the current has the expression ( |3.5| ), this implies a = and that 
non-zero values for a give finite-size corrections. Using Eqs. (|6.5|) - (|6.8|) we get: 



a: 



s? -£ riJ (2(A + l- S ) + H(„-l-2A) 



(6.9) 
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We define 

B= l 1 r g + " ) . (6.10) 

2(l-g + A) v ; 

As will be seen, B will represent the large L limit of the length b of the fluid. Denoting 

v v , . 

7= 5 -j (6.1D 

and using Eq.(|6.9|) we obtain 



., 2(1 + A -q)_J, l-q + 2X 



01 = i* 1 " on ■ ^ 7 • (6-12) 

9-1 V 2(1 — g + A) / 

Eq. ( |6.12|) connects the value of the current and the width of the fluid. We are going to 
obtain a second equation relating the same quantities examining the condensate. We 
write the expression of the current in the condensate in the mean-field approximation 
in a symmetric way 

J = \ (p+(y - j)p-°(y) + P+(y)P c °(y + ±)) 

-\ (p c -(y - \)p+{y) + p-(y)p+(y + j)) (6.13) 



and define 



P±=l±w- (6-14) 

Eq.( |6.3|) is then satisfied for any w. Using Eqs. (|6.13| ) and fl6.14|) and going to the 
continuum limit we find the following expression for the current: 

^ (5 _ 1)( i_ mV £ + I- (6 , 5) 

(an analogous expression for the current appeared for the first time in Ref. [|18|]) The 
same relation could have been obtained directly from Eq. ( |B.3| ) We keep in mind that 
charge conservation gives 

[p±&y = \- (6.16) 

We will now impose boundary conditions for the differential equation ( |6.15| ). We will 
require that the concentration of positive particles is continuous at the end of the 
condensate (y = a) and that of the negative particles at the beginning of the condensate 

(y = o): 

p c » = pI 

P C -°(0) = P R (6.17) 

Notice that in this way, the density of positive particles p+{y) is continuous at the point 
a but discontinuous in the origin where one has a shock even at finite L. Conversely, 
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the density of negative particles p-{y) is continuous for y = but not for y = a. 
Eq.( |6.15 ) and ( 6.17 ) we get the profile of the concentrations in the condensate: 



a 



w 



tan 



where a is given by Eq. 

q + l 

The boundary conditions give: 
,f3a q — 1 



and (3 is related to a: 



P tan( 



-L 



Using 



(6.18) 



(6.19) 



(6.20) 



g + 11 

The equation ( |6.20| ) relates again the current to the width of the fluid. In this way, from 
the Eqs.( |6.12| ) and ( |6.20[ ) we can derive both the current J and b as functions of v, A, q 
and L. Once they are known, using ( |6.6| ) and (|6.18 ) one can derive the concentrations 
in the condensate and the fluid. Knowing the concentrations, one can compute any 
correlation function like co,o, c+ ( _, etc. 

Notice that for y = the density p + is discontinuous even for finite L. Therefore 
the bump ends with a shock on the left. As we are going to show shortly the bump will 
end with a shock also on the right in the large L limit. 

In order to solve the transcendental equations ( |6.12j ) and (|6.20| ), we will do a large 
L expansion. It is convenient to define a new variable u instead of f3: 

2u 



P 



71 



and instead of Eq. (|6.20|) we obtain 
tanw 



2u 



5(1 



a)v 



-i 



where 



S 



(6.21) 



(6.22) 



(6.23) 



La(q - 1) 

Keeping only the term of order u in the left-hand side of Eq.( |6.22|) and using Eqs. (|6.19| ) 
and ( |6.21| ) we obtain: 



a 



5 2 (1 



4(1 



VTT 



-S) 



We now use Eq.( p,12|) to obtain: 

5 2 (q-l) 
7 2(1 + A -q) 



(6.24) 



(6.25) 
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Figure 8. The concentrations of positive and negative particles in a shock (q = 1.2, 
A = 1, p = 0.2, L — 400) obtained from Monte-Carlo simulations and from the model 
as described in the text (solid curves). 



Collecting all the results together, we obtain finally: 

4/ 1 , / (q + l)n V / A(q + 1)(1 - A) 



1 - y \ \ . (6.26) 



and 



where 



q-1 \(q-l)ALJ \ L(q-l)vA 



b-B_(B\ 2 (g + l) 2 7r 2 

V \v) 2(l + X-q)(q-l)A 2 L 2 1 ' ' 



A=l-B. (6.28) 



B is the length of the fluid in the thermodynamic limit and A is the length of the 
condensate. From Eq.( |6.26 ) we learn that in the large L limit the current is indeed 



given by Eq. (|3.5|) and that the finite-size correction terms are not only q dependent but 
(through A) are dependent also on the density of vacancies v and on A. Notice that the 
finite-size corrections to the current although quadratic in 1/L (the relevant parameter 
is 5 given by Eq. (|6.23 )) can be large even for large lattices, say L = 200. We remind the 



reader that in Sec|| we have shown that in the grand canonical ensemble the finite-size 
corrections are linear in 1/L. A different finite-size behavior for the two ensembles was 
already noticed for the disordered phase (see Eqs. ( |3.7| -10)). 



6.1.2. Comparison of the mean-field predictions with results obtained by other methods. 
Finite-size scaling. We are now in the position to compare the result of the model 
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with the Monte- Carlo data. We first start with the current (see Fig.|j). We notice 
that for large values of L {L > 150) the mean-field values obtained using the method of 
Appendix B which are already almost identical with those obtained in the continuum 
approximation given above are close to the Monte-Carlo data. One can be more precise: 
Using Eq. flO&D with q = 1.2 A = lp = 0.2 one gets J = 0.0534 for L = 200 and 
J = 0.0522 for L = 400 to be compared with the values J = 0.059 respectively 
J = 0.053. This observation is interesting since one wouldn't expect a mean-field 
calculation to control so well finite-size corrections. For smaller values of L (see again 
Figl4|), the Monte-Carlo data (confirmed by the calculations done using the algebraic 
approach) as well as the mean-field values show a strange behavior of the current in 
the region 50 < L < 70. The explanation of this phenomenon is simple. For given 
values of the parameters q, A and p for which in the thermodynamical limit one is in the 
mixed phase, the bump picture appears only if the size L of the system is larger than a 
certain minimal value. We would like to add that we have not checked for which values 
of L the continuum version of mean-field coincides with the discrete version presented 
in Appendix B. 

The phase diagram is at this point known. The mixed phase, for a given values of 
A and p, is between 1 < q < q c , where q c is 

*=l + ^. (6.29) 
This value is obtained from the condition 

B = 1 (6.30) 
which says that one has no condensate. 
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Table 1. The estimates of the critical value of q obtained using the grand canonical 
ensemble (gf c ) or Monte-Carlo simulations (q£) compared with the values q c given by 



formula ( gjg ). 



A 


p = m 


= 0.1 


p = m 


= 0.2 


P = 


A A 

m = 0.4 


qf 


q<- 


qf 


Qc 


qf 


q c 


0.25 


1.10 


1.08 


1.2 


1.14 


1.3 


1.22 


0.50 


1.15 


1.17 


1.3 


1.29 


1.6 


1.44 


0.75 


1.25 


1.25 


1.5 


1.42 


1.8 


1.67 


1.00 


1.35 


1.33 


1.62 


1.57 


2.2 


1.89 


2.00 


1.60 


1.67 


2.2 


2.14 


2.9 


2.78 


3.00 


1.85 


2.00 


2.5 


2.71 


3.5 


3.67 


4.00 


2.0 


2.33 


2.8 


3.28 


3.8 


4.56 



One can ask to which extent the mean-field prediction (|6.29|) is correct. In Table [3] 
we present the estimates of the critical points q 9 c c obtained from the "discontinuities" of 
the derivatives of the current and of the two-point functions c_ j+ (l) as a function of q 
using the grand canonical ensemble (see SecfJ) together with the values of q c obtained 
using formula (|6.29|) . The estimates have errors which are hard to evaluate, therefore 
an optimistic point of view would suggest that ( |6.29| ) is exact. A pessimist would say 



that mean-field is a good approximation in the vicinity of the phase transition only for 
low values of A and densities. An argument in favor of the optimistic point of view will 
be given below. 

A well known question in equilibrium statistical physics is the problem of finite-size 
scaling corrections in the determination of critical points JT^|. One can ask a similar 
question here. For a finite value of L (at fixed A and p) an estimate of q c can be 
obtained for example using the fact that in the thermodynamical limit the derivative of 
the current has a a jump at q c which at finite values of L gives an obvious way to get 
an estimate for q c . We would like to know how the estimate approaches q c . In order 
to get an insight to this problem, one can assume that the current changes when the 
condensate appears. Therefore let us go back to Eq. ( |6.271 ) which gives the length of the 
condensate for finite L. One sees that b > B. In particular, for a given value of L, there 
is a value of q that we denote by q c {L) for which there is no condensate, i.e. 

b(L) = l. (6.31) 

Denoting 



q c (L) =q c ~t 



(6.32) 
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and using Eqs. fl6.27| ) and ( |6.31| ) we obtain: 



2 



(2-^(1 -.) (2 -" + A(1 -"" 2 (6 ' 33) 

From Eq.( |6.33| ) we see that e > 0, which means that for a given value of q < q c up to 
a certain length L, one has no bump, in agreement with the interpretation given above 
to the behavior of the current as a function of L in Fig.^. On the other hand Eq. ( |6.33| ) 
gives obviously a finite-size scaling critical exponent 2/3. Unfortunately, this derivation 
is only heuristic. Here are the reasons. A careful examination of the Eqs.( p.l8| - |6.2Cl| ) 



shows that the value b(L) = 1 is never obtained although the large L approximation 
( |6.271 ) gives it. Mean-field theory suggests that the condensation appears directly with 
a non-zero value of the length a and implicitly a value of b < 1. Using this value of b 
(which is obtained from a transcendental equation) one can estimate numerically the 
critical exponent. It is also 2/3 but the right hand side of Eq. (|6.33|) is different. 

A careful examination of the Monte-Carlo data shows a similar phenomenon 
although the value of a at which the condensate appears is slightly different than the 
one predicted by mean-field. 

Instead of working at fixed v and A, one can take fixed q, A (in the mixed phase) 
and L and look for at which density v c (L) starts the condensation (b(L) = 1). We use 
again Eq. (|6.27|) and get 

V v c ) (q-l)(l-q + 2X) y 

where v c is obtained from Eq. ( |6. 1Q|) with B — 1: 

v c = f- q+ J? (6.35) 
l-g + 2A ( J 

From Eq. ( |6.34| ) we learn that for finite L the condensation occurs at a lower density 
of vacancies (larger density of particles) than at the critical point. This behavior can 
also be seen in Fig.[7| where the p dependence of the current at fixed q = 1.2 and 
A = 1 is shown for various lattice sizes. For large values of p one is in the mixed phase 
( J = (q— l)/4). If we decrease the value of p the current increases (see Eq. (|6.26| )) up to 
a point where the condensate disappears. For small values of the density the mean-field 
approximation does not apply (one is in the disordered phase). If one increases the size 
of the lattice, the minimal value of p where mean-field still applies gets smaller and the 
current gets closer to its asymptotic value (g — l)/4. One way to estimate v c is to take 
the values of the density p c where the current has a maximum. One can then estimate 
the exponent alpha from the expression 

1 - = const.L-* (6.36) 



33 



Using data taken at L = 100, 200, 400, 800 and 1600 one gets a = 0.5 ±0.2 which makes 
the exponent a = 2/3 given by Eq. (|6.34 ) plausible. Moreover, this also implies that the 
mean-field value for v c given by Eq. ( |6.35| ) is probably correct. The estimate of a has 
large errors since it is difficult to determine precisely the position of the maximum of 
the current. 

We are going to verify the mean-field predictions in the mixed phase away from q c . 
We take again q = 1.2, A = 1, p = 0.2 and L = 400. Solving Eqs.(|BTT2D and flggg ) 
numerically and using Eqs.( |6.l8D and Q6.19| ) we get: 

J = 0.05257, a = 0.3200, w = 0.1134 tan ^8.244(^ - y)\ (6.37) 

We have done Monte-Carlo simulations taking only configurations in which a had 
approximately the value given in (|6.37 ). For these configurations we have determined 
the concentrations of positive and negative particles. They are shown in Fig.[| together 
with the prediction obtained from Eq. fl6.37|) . The agreement couldn't be better. One 
can use (|6.37| ) in order to compute the correlation functions co,o and c+ _ and compare 
them with the Monte-Carlo data. This comparison can be seen in Fig.|6| and was repeated 
for all correlation functions. The agreement is as good as for cq+ and c_ + . 



6.1.3. Profiles for large L We would like to discuss now the composition of the bump 
in the thermodynamical limit. We denote the concentrations by R instead of p. From 
Eqs.flOg), (fT24j) and (|6T?D , we get 

= 5, < = p Rl= 1 -^ (6-38) 

This result is mostly interesting. One can see the closed system as two open 
systems, one containing only positive and negative particles (the condensate) with 
equal concentrations and a second one containing also vacancies (the fluid). From 
our knowledge of the two-state model |J with the bulk rates = q, = 1 we 
know that the current has the expression (|3.5|) from mean- field (it is just a repetition of 
the calculation done above where Eq. (|6.17|) is replaced by an equation which takes into 
account the boundary rates). This is the so-called maximum-current phase || [2|. In 
the same maximum-current phase but only for the case of totally asymmetric diffusion 
it was shown that the correlation functions have an algebraic decay. It is plausible to 
assume that the same applies for the partially asymmetric case which is ours. This 
would also suggest that the fluid phase is also algebraic. This assumption can be 
checked independently since the algebraic tools for this problem are already known 
P0|| . Unfortunately we were not able to check if indeed the correlation functions c 0j o, 
c + - etc. in the present model have an algebraic decay. This would imply to work 
not with the macroscopic variable y but with R (see Eq.( [4.3| )) and to take lattice sizes 
beyond our means. 
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The reader might have noticed that all the expansions used above diverge for q — 1 
[see Eq.( |6.23| )). We proceed now to discuss this case. 



6.2. q = l 

A repetition of the calculation done for the case q = 1 gives the following results. The 
current in the fluid phase has the expression: 

A . v , v , 

2 ~ ih < 6 ' 39 > 

In the condensate, the current is 
1 dw 

J = -I*7 (6 ' 40) 



which can be obtained from Eq. fl6.15|) for q — 1. Taking into account the boundary 
conditions, we obtain: 



and 



«-=£<!-„) (6.41) 



Equating the two expressions for the current (|6.39| ) and (|6.42| ), we find b as a function 



of L and consequently, J as a function of L. In order to show how good our results are, 
one can compare for A = 1, p = 0.2 and L = 100 (which is a small lattice) the value 
of the current obtained using ( |6.39|) - (|6.42|) which is J = 0.02588 with the Monte-Carlo 



value J = 0.0256. 

In Fig.||, we compare the profiles, now for q — 1, just as we did in Fig.|8| for q > 1. 
Notice the nice straight lines. 

For large values of L one obtains the following leading contributions to the current 

and to the length of the fluid 

aZJl -v) 



V \ 1 + TT7V \) ( 6 - 44 ) 



which explains Eq.( [3.10| ). 

In the thermodynamic limit, keeping the notations used already, we get: 

B = v, A — l — v (6.45) 

Rl = 0, B% = l-!L (6.46) 

i? fl = o, R°° = ?L (6.47) 
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and of course J = 0. 

It should be interesting to study in detail the properties of the fluid and the 
condensate at q = 1 looking not at macroscopic distances (y) but at microscopic ones 

6.3. q < 1 (pure phase) 

One can repeat the calculations done for the mixed phase for the pure phase. We give 
directly the results. Instead of Eq.( |6.18| ) one has: 



V2 /„ i / at \ V2 




which gives the profiles of the concentrations in the condensate, the boundary condition 
( |6.2(]| ) been replaced by 



w (0) = pi = v/b (6.49) 
This allows to obtain for large values of L the expression of the current: 

X+l — q \q + l J 
and of the density of vacancies in the fluid: 

pi = 1 - 2 J/A (6.51) 

the last expression shows, as expected, that in the large L limit, the densities of vacancies 
in the fluid tends exponentially to 1 which is precisely what one expects in the fluid 
phase. The expression of the current given by Eq. (|6.50p is compatible with the fits done 
to the Monte-Carlo data (see for example ( |3.9| )). An exponential fall-off of the current 
was also observed in a different model [Ell. 



7. The spontaneous breaking of translational invariance. Flip times in the 
pure and mixed phase 

In the last sections we have described the pure and the mixed phases. Both phases show 
long range correlations which are different from a quantitative point of view (exponential 
versus algebraic finite-size convergence etc.), but up to this point we haven't proven that 
they are also qualitatively different. Moreover, mean-field calculations (see [Appendix B| 



and Sec.0) show no qualitative differences between the two phases: one obtains pinned 
configurations in both cases. We are going to show now that half of these mean-field 
predictions are wrong and that in the pure phase translational invariance is broken and 
that this is not the case for the mixed phase. 
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Figure 10. The Ginzburg-Landau potential in the pure phase (p = 0.2, A = 1, 
q = 0.9) as a function of \M\ obtained from Monte-Carlo simulations for L = 100(o), 
200(+) and 400(d). 
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Figure 11. The Ginzburg-Landau potential in the mixed phase (p = 0.2, 
A = 1, q = 1.2) as a function of \M\ obtained from Monte-Carlo simulations for 
L = 200, 400, 600, 800, 1000 (from top to bottom). 



It is useful to consider the order parameter M denned as follows: 
k=i 



(7.1) 



One notices that the order parameter is sensitive only to the vacancies and to the 
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Figure 12. The Ginzburg-Landau potential in the disordered phase (p = 0.2, 
A = 1, q = 2.5) as a function of \M\ obtained from Monte-Carlo simulations for 
L = 400, 600, 800, 1000, 1200, 1400 (from top to bottom). 



macroscopic properties of the system (it sees only the lowest frequency). For a given 
length L of the system we have measured using Monte-Carlo simulations the probability 
P(|M|) to find a certain value of the modulus of M. Next we have considered the 
Ginzburg-Landau potential 

f(\M\) = ~logP(\M\) (7.2) 

In Ref.0] we have proven, for another order parameter and another model, that this 
function converges in the thermodynamical limit. In Figs.|TU|, |TT| and |T2| this function is 
shown for the pure, mixed and disordered phases. 

These figures confirm know results. Indeed if one has a block of vacancies of length 
b (again macroscopic length!) a trivial calculation shows that f(\M\) should have a 
minimum at the value 

|M| min = 4sin(7r6) . (7.3) 
no 

Here v is as usual the density of vacancies. In the pure phase (b = v) this gives the 
value |M| min = 0.3027, for the mixed phase one can use the value b = 0.675 obtained 
using Eq. ( |6.10| ) which gives |M| min = 0.238 and, finally, for the disordered phase one 
gets |M| min = since here b = 1 in agreement with the values seen in the figures. 

We would like at this point to make, what is in our view, an important observation. 
In Ref . |J it was advocated that the Ginzburg-Landau potential defined in Eq. ( |7.2| ) is a 
useful concept not only in equilibrium statistical physics but also in the present context. 
As one can see from Figs.[TT| and this is certainly the case for the mixed and disordered 
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Figure 13. The average flip time T defined in the text as a function of L in the pure 
phase (p = 0.2, A = 1, q = 0.8) obtained from Monte-Carlo simulations. 
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Figure 14. Same as in Fig.[l3| for the mixed phase (q — 1.2). 



phases since this function nicely converges for large L. This is not the case for the pure 
phase. Notice that in this case the function which converges is rather f(\M\)/L. 

A similar phenomenon occurs for another concept "exported" from equilibrium 
statistical physics to non-equilibrium, namely the position of the roots of the partition 
function seen as a function of fugacity p2] . As shown in Ref. p3| , they converge to a 
parabola in the disordered phase, to an ellipse in the mixed phase, but shrink for large 
L for the pure phase. These observations might be relevant for other models or for the 
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same model in two dimensions since it might give a way to distinguish between the pure 
and mixed phases. 

In order to obtain new information one can look not to \M\ but to M itself and, 
instead of considering the stationary state only, examine the time dependent behavior 
of the system. Using Monte-Carlo simulations we have looked for a given size L of 
the system, at the average flip time T from a configuration which is in the domain 
Re(M) > 0, -0.1 < Im(M) < 0.1 to a configuration Re(M) < 0, -0.1 < Im(M) < 0.1. 
Here Re(M) and Im(M) represent the real respectively the imaginary parts of M. These 
average flip times are shown as functions of the size of the system L in Fig.|l3| (for the 
pure phase) and in Fig.|l4] (for the mixed phase). For the pure phase we limited ourselves 
to L = 140 since for larger values of L the flip times are to long. Now we notice a major 
difference between the two phases: whereas in the pure phase the flip time increases 
exponentially with L (exponentially increasing barrier), in the mixed phase the flip time 
increases algebraically. Actually the data are compatible with a linear dependence. This 
result implies that in the pure phase but not in the mixed phase, translational invariance 
is spontaneously broken. This settles the problem about the difference between the pure 
and the mixed phase. 

If one uses as a definition of charge segregation the existence of correlation functions 
which are dependent of the macroscopical variables (y in our notations for the two-point 
functions), we have to distinguish between the two possibilities put in evidence in this 
section. One in which one has the breaking of translational invariance, which we propose 
to call "charge segregation of type B" and one in which translational invariance is not 
broken which we propose to call "charge segregation of type UB" . 

We would like to mention that another approach to average flip times in order to 
put in evidence charge segregation can lead to ambiguous results. A different model was 
studied in Ref. |24]]. The authors looked, (this would correspond to our model), at the 
average flip time T from a configuration with &i|M| min to another with ^l^lmm- Here 
b\ and 62 are L independent constants. They found that T increases exponentially with 
L. From the shape of the Ginzburg-Landau potential shown in Figs.|l0] and |Tl| the same 
should be true also in our case for both the pure and mixed phases. We could distinguish 
between the pure and mixed phases only by considering the average flip times in the 
complex M plane. 

We can now continue the analysis of the model and consider the disordered phase. 



8. The disordered phase. 



For fixed A and densities and q > q c where q c is given by Eq.( |6.29|) , one finds the 
disordered phase. In this phase the correlation length is finite. What is the information 
we have already about this phase? 
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Figure 15. The connected correlation function c+—(k) as a function of k for for 
p = 0.2, A = 1, L = 100 and three values of q. The data are computed using the grand 
canonical ensemble. 
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Figure 16. The inverse correlation length as a function of q for p — 0.2, A = 1. The 
data are obtained determining the correlation lengths from systems of various lengths 
L (up to 200) and then extrapolating the values to their thermodynamical limit. The 
dashed curve is given by Eq. 



a) For large values of q, the current is known (see the large q expansion presented in 
Sec. 2) and is given in the case of equal densities of the charged particles by Eqs.( |2.47|) , 
( gg§ ) and ( P3BQ . 

b) For A and q related by the relations fl2.i2| ) (those are curves in the g-A plane 
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Figure 17. The density p as a function of the fugacity z for q = 1.2, A = 1, 
L = 100, 150, 200 sites. 



in the domain where one is in the disordered phase), the quadratic algebra (|2.3| ) has 
finite-dimensional representations. If the representation is finite-dimensional, it follows 
that the correlation length is finite. The calculation of various correlation functions is 
simple but nevertheless tedious. We have not done any calculation of this kind because 
we didn't expect any surprises. 

Since we are interested to understand the nature of the phase transition at q c , we 
have taken A = 1, p = 0.2 (for this case we know, using Eq.( |6.29 ), that q c = 1.57) 



and computed the correlation function C_| clS db function of the microscopic distance k 

for three values of q. We have taken L = 100 and used the grand canonical ensemble. 
The results are shown in Fig.|l5| in a double-logarithmic scale. One sees that in the 
disordered phase (q = 1.7 and 1.8) one has an exponential fall-off but that at q c one has 
an algebraic one. A fit to the data gives: 

c+ _(jfe) = (0.04 ± 0.01)£r L2±a3 (8.1) 

Next, we were interested in the q dependence of the correlation length. For this purpose, 
for a given value of q, we have, using the grand canonical ensemble, taken lattices of 
different lengths L, and determined the correlation length £ fitting £ of equation 

c+,_(Jfe) = A AT 1 exp(-k/£) (8.2) 

The factor k^ 1 was chosen in order to improve the fits to the data. Then we have 
extrapolated the values of £ to obtain the large L limit of £. These values are shown 
in Fig.|TBl As expected the inverse correlation length vanishes at the critical point. As 
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seen from the figure, the data are compatible with the fit: 

r 1 = 0.64(g -g c ) ' 6 (8.3) 
At q c , the correlation length has the following large L behavior: 

f = 0.0065L + 3.98 (8.4) 

Taking into account that one obtains an excellent fit to the data, the readers who like 
conformal invariance should give a second thought to this result. 

It is obvious that more work is needed in order to fix the values of the two exponents 
in Eqs. (|8.1| ) and (|8.3| ), check for universality, etc. Our attention however was taken by 
a problem which we have considered more interesting, namely to understand which 
symmetry is spontaneously broken at the phase transition. We remind the reader, that 
the symmetries of the system are translational invariance (for which we have shown in 
Sec.|7| that it is not broken in both the disordered and mixed phase), CP and U(l) x U(l) 
corresponding to the conservation of the numbers of positive and negative particles. The 
answer to this question is given in the next figure. 

In Fig.|l7] we show for A = 1 and q = 1.2 (which for p = 0.2 corresponds to the 
mixed phase but becomes q c if p = 0.055), the fugacity(,z) dependence on the density(p) 
for three values of L. One notices that the L dependence of the curves suggest that in 
the thermodynamical limit, p increases linearly with z up to the point where p reaches 
the value 0.055 (the linear increase is an empirical observation). At the corresponding 
value of z, one gets a vertical and z does not fix anymore the density. Such a behavior 
of the density is known in Bose-Einstein condensation |25j where (there one has a U(l) 



symmetry) the conservation of the number of particles is broken. In the next section we 
are going to discuss this topic in detail. 



9. Spatial condensation 

As is well known, Bose-Einstein condensation is a quantum mechanical phenomenon in 
which identical particles condensate in a state of zero momentum. This phenomenon 
takes place in three or more dimensions. In a simple but fascinating one-dimensional 
stochastic model with a "defect" (reviewed in Appendix C), M. Evans |26j has shown 



that, choosing the right variables, the probability distribution for the stationary state 
can be written in terms of occupation numbers with a "dispersion" relation (they are 
no momenta involved) which give the equivalent of a Bose-Einstein condensation (there 
is a zero bosonic mode). In this model the condensate corresponds to an empty domain 
(no particles, just vacancies). 

We will pursue the investigation of our model in order to see if indeed more properties 
are common to the usual Bose-Einstein condensation and to the transition at q c . The 
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Figure 18. £ defined by Eq.(|J) as a function of L for p = 0.2, A = 1 for q = 0.9, 1.4, 
1.6 and 2.2 (from top to bottom). The data are obtained using the grand canonical 
ensemble. 
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Figure 19. The "pressure" P defined by Eq.({3.3|) as a function of 1/p for q = 1.2 
and A = 1. The data are obtained as explained in the text. 



compressibility k is defined by the relation: 

k = ( 2 L 

where 

<(A /3 ) 2 > 1 /2 



c 



p 



(9.1) 



(9.2) 
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Figure 20. The Ginzburg-Landau potential /(«;) as a function of the local density of 
vacancies vi. The curves are for q = 1.6 « q c and L — 400(a), L — 800(6), L = 1600(c), 
and for q = 2.2 and L = 400(d), L = 800(e), L = 1600(/). In all cases we have taken 
A = 1 and p = 0.2. 



Using Eq. (|2.36|) we have measured ( as a function of L for A = 1, p = 0.2 and several 
values of q. The data are shown in Fig.|l8[ We notice that for q = 0.9 (in the pure 
phase) and q = 1.4 (in the mixed phase), ( tends to a constant which gives a divergent 
compressibility. This is not the case in the disordered phase (q = 2.2). At the critical 
point (q c = 1.6) ( does not decrease like L~ 1 / 2 but the lattice sizes are too small in order 
to decide its ultimate behavior. 

Another quantity of interest is the "pressure" P, which can be defined as for 
equilibrium statistical physics: 

P=ilogZr 

and can be determined using Eq. ( |2.41| ) 



(9.3) 



J = Xz 



Z L-1 



Az(l-— logZn 



(9.4) 



Assuming that log is an extensive quantity one obtains: 

In order to find an "isotherm" (there is no temperature in our problem) we have, for 
q = 1.2 and A = 1, computed the current as a function of p (at each value of p we 



(9.5) 
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have have computed the current for various lattice sizes and extrapolated to large L) 
The "pressure" was computed using Eq. 



and is shown in Fig. 19. There are several 
interesting things that we can learn from this figure. First that the "pressure" is positive, 
next that it decreases when the density gets smaller. There are no theoretical reasons 
that we can think of that can justify why the "pressure" defined by Eq.(IO) should have 
the same properties as the pressure in equilibrium systems. Moreover, the "isotherm" 
shown in Fig.lIU] looks similar to the isotherm of the free Bose gas when Bose-Einstein 
condensation takes place. 
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Figure 21. The correlation functions c +i _ for the canonical ensemble (bottom, 
computed by Monte-Carlo simulations) and for the grand canonical ensemble (top) in 
the mixed phase: A = 1, q = 1.2, p = 0.2, L = 100 and L — 200 sites. 

Having in view other systems where spatial condensation might take place we have 
looked for another method to put it into evidence. We have looked at a given lattice 
size L at the probability to find a local value of the density of vacancies V\ which is 
defined as a local average in the macroscopic y coordinate (vi = number of vacancies in 
the interval Ay divided by Ay). The corresponding Ginzburg-Landau potential f(v{) 
for various values of L with q in the disordered phase and at the critical point is shown 
in in Fig.[2(| In the disordered phase (q = 2.2) f(vi) has, as it should, a minimum at 
v = 0.6 which corresponds to the average density. For q = q c the function f(v{) is almost 
flat which would give a Ginzburg-Landau potential corresponding to a first-order phase 
transition although we know that the transition is second-order. A similar situation 
occurs also in the case of Bose-Einstein condensation where some thermodynamical 



quantities behave like in a first-order phase transition [25 



At this point one can ask oneself to which extend the grand canonical ensemble is 
useful at all in the mixed phase (we have used it consistently to compute the current). 
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Figure 22. ( defined by Eq. (|9.2| ) as a function of L for p — 0.2, A = 1 and q = 1.2 
for various values of v (from top to bottom v = 0, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1) 
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Figure 23. The correlation function c+,_ as a function of y for p — 0.2, A = 1 q = 1.2, 
L = 100 and various values of v = 0,0.05,0.06,0.07,0.08 (from top to bottom). The 
solid line corresponds to the fluid (see the text). 



It turns out that if one wants to compute local quantities like the current or other local 
two-point functions (see Fig.[|), the grand canonical ensemble gives correct results. We 
have to keep in mind that the finite-size effects are of order (see Eq.( |3.9| )) as opposed 
to L~ 2 for the canonical ensemble (see Eq.( |6.26 )). If, however, one is interested in the 
values of the correlation functions at fixed macroscopic distances y (see the definition 
in Eq.(f4.3p) using the grand canonical ensemble one gets wrong results. In order to 
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illustrate this point, in Fig.^l] we compare the values of the correlation function c +i _ 
obtained using the grand canonical ensemble and the canonical ensemble (Monte-Carlo 
simulations) in the mixed phase. We notice that around y = or 1, the two methods 
give the same results but that the results differ drastically otherwise. 

Can the grand canonical still be used in order to get information about the mixed 
phase? We believe that the answer is yes and we have explored this possibility. The idea 
is inspired by the way the problem is solved for the usual Bose-Einstein condensation (the 
procedure is shortly repeated in Appendix C). There, one introduces in the Hamiltonian 
a supplementary term, which breaks the U(l) symmetry of the problem, proportional 
with a constant v and the square of the volume. After a Bogoliubov transformation, 
the net effect is a modification of the partition function given by the grand canonical 
ensemble by a factor which depends on the volume, v and the chemical potential. This 
changes the relation between the density and the chemical potential. In order to compute 
a physical quantity, one keeps v fixed, then one takes the volume to infinity and at the 
end, makes v equal to zero. Since in our case we were not able to find a way to break 
explicitely the U(l) x £7(1) symmetry of the problem, we have, in analogy with what 
is done for the Bose-Einstein gas, just modified the the partition function ( |2.32| ) of 
the grand canonical ensemble by a factor which takes into account that we have two 
chemical potentials and that at the critical point the chemical potential is not zero: 

Z^\v) = exp(- — )Zj?>» (9.6) 

2fX c - fix - fl 2 

In Fig.^2| we show the L dependence of ( defined by Eq. (|9.2|) in the mixed phase for 
various values of v. We have to keep in mind that the range of values of L is limited. For 
v = 0, 0.05, 0.06 and 0.07 it looks like the compressibility k (see Eq. (|9.1| ) diverges. A 
dramatic change takes place for v = 0.08, the compressibility becomes convergent. It is 
plausible to assume that the compressibility converges for any non- vanishing v (one has 
to take larger lattices to observe the phenomenon). We further investigate the effect of 
the modification of the partition function. In Fig.^3| we show again the the correlation 
function c +i _ as a function of y (compare with Fig.^T]), but this time for various values 
of v. 

We notice that between v = 0.07 and v = 0.08 again a dramatic change takes place: 
the correlation function (it is not the connected one!) drops to a constant value. It 
is tempting to assume that this should be the correlation in the fluid observed in the 
mixed phase (the symmetry breaking term eliminates the condensate also in the case of 
the Bose-Einstein gas). This would imply: 

c+ _ = Rl = (i?A )2 (9.7) 

where is given by Eq.( |6.38D . Using Eq. ( |6.10| ) this gives c +> _ = 0.031. This is the 
solid line shown in Fig.BS] which gives approximately the correlation function observed 
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for v = 0.08. 

It is obvious that this investigation is only preliminary but we believe that it gives 
the right approach to derive the properties of the fluid in the mixed phase. It also 
shows how similar are the condensation phenomenon in the Bose-Einstein gas and in 
our model. 



10. Conclusions 

We believe that we brought some more light in the understanding of phase transitions 
in one-dimensional stationary states. The model we considered has two advantages: it 
has interesting physical properties and has behind it a mathematical structure given by 
a quadratic algebra. This mathematical structure is not only a useful computational 
tool but allows, due to its similarity to the mathematical formulation of equilibrium 
statistical physics, to define relevant quantities which characterise the nature of the 
phase transition. A first example is the definition of the grand canonical ensemble. 
Working with this ensemble has not only allowed to obtain exact results (like the 
large q expansion presented in Sec. |2.4|) but also to do precise numerical calculations 
on large lattices. Moreover, it allowed to define the analogue of the pressure which, 
amazingly, has the same properties as in equilibrium (see Fig.[T9|) . Using the grand 
canonical ensemble we were able to identify the nature of the phase transition between 
the disordered and the mixed phase since this transition resembles the one seen in the 
Bose-Einstein condensation (see the fugacity dependence of the density in Fig.[T7] and the 



volume dependence of the compressibility in Fig.|T8|) . Most importantly, using the grand 
canonical ensemble, one can break explicitly the continuous symmetries of the problem 
(conservation of the number of charged particles) in analogy with the treatment of the 



same problem in the case of Bose-Einstein condensation (see Eq. fl9.6|) and Figs.|22] and 



Before discussing the physics of the model, we would like to stress some other results 
presented in this paper. We have shown that the Ginzburg-Landau potential, introduced 
earlier in Ref. || to study first-order phase transitions, is a meaningful quantity also in 



the present context. This was checked for two order parameters (see Figs. 11, 12 and 



20]) and the behavior of this function in the case of a second-order phase transition was 
for the first time observed. The Ginzburg-Landau potential can also be useful in order 
to put in evidence the breaking of translational invariance since in this case it does 
not converge in the infinite volume limit. This observation is based only on the data 



presented in Fig .[10] and obviously one needs to do more work to clarify this point. 

We think that we have also found the proper method to look at flip times in order 
to distinguish between situations when one has charge segregation with and without 
translational symmetry breaking. As described in Sec.[7|, one has to look not only at the 
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module of the order parameter defined by Eq. ( |7.1[ ) but also at its phase. 

Inspired by the numerical results presented in Appendix B , in Sec.|| we have shown 
in detail how to do mean-field calculations for this model. The present method can 
also be used for a larger class of models. These mean-field equations are two coupled 
non-linear differential equations (see Appendix B ). The fact that one is able to find 
analytic solutions for the stationary states is very non-trivial. We will return to this 
topic in Ref . |7| where we consider the case of different densities of positive and negative 
particles. We will explain why the differential equations should be called the two- 
component Burgers equations and give more analytic solutions of these equations. 

We found the physics of this simple model fascinating. There are two phase 
transitions. One separates the disordered and mixed phases and is of second-order. 
This phase transition is related to spatial condensation. We got estimates for several 
critical exponents related to the inverse correlation length, the algebraic behavior of the 
two-point function at the phase transition and finite-size scaling corrections (see Secs.^| 
and [8p. The mixed phase with its simple expression for the current and its independence 
on the density and other parameters as well as the bumps well described by mean-field 
are nice pieces of physics. The nature of the bump is also interesting: one has two 
open systems glued together. One, the condensate composed of charged particles is well 
understood from the corresponding two-state problem. About the second, the fluid, we 
know little except for the fact that it shows no macroscopic structure. Actually, one 
major difference between Bose-Einstein condensation and spatial condensation is that 
here the condensate has a simple structure whereas the fluid phase is interesting since 
it is different than the disordered phase and hence worth studying. 

In another publication |E7| where we consider the open system with the same bulk 
rates as in the present one we will show that a first-order phase transition occurs between 
one phase which is the condensate and a second one which is the fluid. In this way one 
can get directly more information about both phases. 

Unfortunately we did not pay much attention to the second phase transition, the 
one between the mixed and pure phase, where translational invariance is spontaneously 
broken. 

This work has left us with quite a few open questions which we would like to address 
in some future publications: 

(a) We have not touched the study of the model in the case of different densities. A 
second paper will be devoted to this. 

(b) The algebraic approach and the Monte-Carlo simulations were used only for 
a limited numbers of points in the q-\ plane. Most of our conclusions including the 
transition separating the mixed phase form the disordered phase are mean-field results 
which were therefore checked only in a limited region of the q-\ plane. This does not 
mean that in other regions mean-field is valid. 
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(c) The dynamics of the model with the parameters chosen such that the stationary 
states correspond to the pure and mixed phases as well as at the phase transition points 
can lead to plenty of interesting results. 

(d) Last but not least, there is many information given in this paper which deserve 
a more accurate research. We did our best but we feel that this is not enough. 
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Appendix A. The connection between the representations of the quadratic 
algebra and of the quantum algebra U q so(2, 1). Recurrence relations. 

The expressions (P~7|) - (|^TTD| ) of the representations of the quadratic algebra ( |2.3| ) suggest 
to any reader familiar with the representations of quantum algebras a connection 
between these very different mathematical topics. This connection is for the purposes 
of this paper not only a mathematical curiosity but, as will be seen, one can write 
recurrence relations for physical relevant quantities which have an interpretation from 
the quantum algebras point of view. We will comment about these recurrence relations 
at the end of this Appendix. 

We first remind the reader of the definition of the quantum algebra U s o(2, 1). This 
is the deformation given by a parameter s of the well known algebra so(2, 1). (The 
deformation parameter is usually denoted by q but this symbol is already used in this 
paper). The quantum algebra is defined by three generators k, e and / satisfying the 
relations |E8|: 



[e,f] 



k 2 - k~ 2 
s — 



kek 1 = s 1 e 

kfk~ l = sf (A.l) 
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One representation of this algebra has the following matrix elements: 



where 
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(A.3) 



(A.4) 



Notice that if 



1 

A 



(A.5) 



then G\ and G2 (see Eq. (|2^ )) can be expressed in terms of fc, e and / as follows: 



k 
s 



e + 



k-k 



-1' 



Go 



k -(f + k - 

s \ s 



(A.6) 



In order to show how the connection between the representation of the quadratic algebra 
and the quantum group can be used, we will consider a simpler example where s = 1 
(this is the familiar case of the so(2, 1) algebra). In this case, one has: 



G 1 = S Z + S~ 



G 2 = S Z + S H 



(A.7) 



with: 



k \k> 



k(k - 1) \k- 1> 



S z \k> 
s- \k> 

S+ \k> = ^k{k + l) \k + l> 
where k > 1. We have to keep in mind (see Eq.i 



expression: 



G \k>= 4,i \k> 



(A.8) 

that the matrix Go has a simple 

(A.9) 
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The calculation of the current for example, in the grand canonical ensemble, implies the 
calculation of the quantity (see Eqs. ( |2.32 )-( [2.34|) ): 

C N \1> (A.10) 

where 

C = G + (z 1 + z 2 )S z + z 1 S- + z 2 S + (A.ll) 
Writing C N |1> in the \n> basis: 



A' 



C N \l>=J2t ( n N) \n> (A.12) 

n=l 

one can find a recurrence relation for the using ( |A.8|) and ( |A.9|) . It is useful to 
denote: 

^*S V) =«f° (A- 13) 
and we get the following recurrence relations for 2 < n < N — 1: 



J^ U N -{Zl+Z 2 )U N + Z 2 U N _ 1 

^«g5'> = *«g0 (A. 14) 

If in the recurrence relations ( |A.14j) one drops the 1 in the term (1 + zi + z 2 ) one obtains 
recurrence relations defined only by the algebra so(2, 1). Similar recurrence relations 
can be found also in the deformed case (s different from 1). We just did not have the 
time to exploit the recurrence relations just mentioned, they might appear however in 
a further publication. 

Before closing this Appendix let us observe that for 

g=- = A + l (A.15) 

r 

G\ and G 2 can be expressed through the generators of the r deformed bosonic creation 
and annihilation operators: 

G x = 1 + VI - r a , G 2 = 1 + Vl-ra + (A.16) 

where 

aa + — ra + a=l (A. 17) 

and 

a |1>=<1| a+ = (A.18) 
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This allows again to write recurrence relations like (|A.14 ) with a clear mathematical 
meaning. 

In defining the grand canonical ensemble (see Sec. 2) we had to make a 
choice between isomorphic algebras (Eqs.( |2T33| )-( |2.34] )) corresponding to different 
normalizations of the generators Go, G\ and G%. As we can see from Eqs. (|A.8| ) and 
( |A.16| ), the chosen normalisation corresponding to Eq.( |2.30| ) looks reasonable. 

Appendix B. Mean-field on the lattice and in the continuum 




Figure Bl. The inhomogeneous solution of the mean-field equations in the pure 
phase for q = 0.9, A = 1, p = m = 0.2 and L = 100. 



Let pk, m k and v k be the average concentrations on the site k of positive, negative 
particles and vacancies. Neglecting correlations, the expression of the currents of positive 
and negative particles on the link between the sites k and k + 1 are: 

J k,k+i = QPkm k+1 - m k p k+1 + Xp k (l - m k+ i - p k+1 ) 

J k+i,k = qPkm k+1 - m k p k+1 + Xm k+ i(l - m k - p k ) (B.l) 

The time evolution of the local densities can be derived from the master equation: 

^ ^ — T+ T+ 
^Pk — J k -l,k J k,k+1 

-77™fc = Jk+i,k ~ J k,k-1 ( B -2) 
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For the stationary state one gets a system of nonlinear equations for the local 
concentrations which might give several solutions. In order to find them, we have 
proceeded as follows. We have used Eq.( |B.2| ) taking discrete values for the time and 
chose the parameters (q, A and p) corresponding to the pure and mixed phase. As 
an initial distribution we have put at random pL positive and pL negative particles 
on the lattice. Starting from one of these configuration the iteration converges either 
to the homogeneous solution (pk = m,k = p) or to inhomogeneous solutions. The 
inhomogeneous solutions have all the same profiles which differ only in their position on 
the ring. Therefore, in mean-field, translational invariance is broken in both the pure 
and mixed phases even for finite L. In Figs. [Bl] and [B2] the inhomogeneous solutions 
are shown for the pure and mixed phase. In both figures one sees two domains. In one 
domain one has only charged particles with non-uniform distributions (the condensate) 
and in the other domain one has uniformly distributed vacancies and charged particles 
(the fluid). For a given values of the parameters q, A and p, the inhomogeneous solution 
appears only for a lattice size L > L min . For larger lattice sizes, one finds inhomogeneous 



solutions different of those shown if Figs.pl] and p2| . Namely if L > 2L m j n one finds 
solutions looking like two pictures shown in Fig |B2| glued together (the same is true 
for the pure phase). Otherwise speaking, for a lattice size 2L one gets a solution with 
two condensates, in each condensate one sees the same concentrations profiles as in the 
case of a single condensate for a lattice size L. The reader can guess what happens for 
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L > 3L min . Based on this observations, we now look at the mean- field equations in the 
continuum. 

Instead of the Eqs.( |B.l|) for the currents, one gets: 
J+ = (q- l)p + p_ + Ap+p 

+ ^7 (?+ 1 )(P+^ ) + A IP+^ P0- 



2L \ dy dy dy dy J 

J- = (q- l)p + p- + Ap_p 

where p + , p_ and po are the densities of positive, negative particles and vacancies, 

P+(y) + P-(y) + Po(v) = 1 (B.4) 
The Eqs.( [B.2| ) are replaced by: 

dp + dJ + dp_ dJ~ 



dt dy ' dt dy 

The solutions of the Eqs.([B.5]) have to be positive, fulfill Eq. (|B.4j) as well as 



(B.5) 



i ,i 

, P+(y)dy = p, / p-(y)dy = m (B.6) 
/o Jo 

Notice that solutions of the mean-field equations will depend on the parameters 

L, q, A, p and m. 

A more detailed study of the solutions of Eq.(|B.5|) will be presented elsewhere M, 



in the present paper we are interested in the stationary solutions which are obtained 
when p = m and with only one condensate. These solutions are obtained in Sec.|| based 
on the observation that in the condensate one has only two species which makes the 
problem solvable and that in the fluid the concentrations are uniform. It turns out that 
in the presence of fluctuations only the solutions with one condensate are relevant for 
the stationary states. 

Appendix C. Bose-Einstein condensation in equilibrium and spatial 
condensation in stationary states. 



We will first review the phenomenon of condensation for the free Bose gas [^9]. The 
partition function Z of the system is: 

Z = Tr exp ((e(fc) - fi)n n - Wu(a + a+)) (C.l) 

k 

where (3 is the inverse of the temperature, e(k) is the energy, k the momentum, p the 
chemical potential, V the volume, v gives the intensity of the symmetry breaking term 
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and is the occupation number: 

n k = a % a k ( C - 2 ) 

The a^'s are Bosonic annihilation operators. We assume for simplicity that we consider 
the three-dimensional gas with a quadratic dispersion relation so that one has Bose- 
Einstein condensation. Since e(0) = we have a zero Bosonic mode. In the absence 
of the symmetry breaking term (y = 0), the Hamiltonian is invariant under U(l) 
transformations 

a ^a e i(t> (C.3) 
where is a phase. Through the non-unitary Bogoliubov transformation 

b = a + W- (C.4) 

the partition function becomes 

Z = e~ — Tr exp -(3 ]T (e(Ar) - fi)n % ) (C.5) 

k 

where are again number operators (nr = 0, 1, 2, . . .) 

% = Hh ( c - 6 ) 

The net effect of the Bogoliubov transformation is the appearance of the chemical 
potential not only under the trace operation but also in the overall factor. This implies 
that p defined by the expression 

"=wi^ z < C7 > 

is the average density of particles only for v = 0. The main point is that we first do 
the calculation for v different from zero, take the infinite volume limit and only at the 
end let v vanish. In this way one gets the thermodynamical properties of the gas also 
under the critical temperature since the condensate (which is a state of zero momentum) 
disappears from the calculation as a consequence of the symmetry breaking term. This 
method can also be seen as a regularization procedure. 

That the free Bose gas is a quantum system (the expression ( p.5|) contains only 
occupation numbers which are c-numbers) can be seen once we compute measurable 
correlation functions like 

g{n,r 2 ) =<^(f 1 )^ t (f 2 )> (C.8) 

where 

m = ^E^*k (c.9) 
vv k 
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since the creation and annihilation operators verify the Heisenberg algebra. That much 
on the Bose-Einstein gas. 

In the case of the spatial condensation which appears in the mixed phase, the 
condensate is separated in space from the fluid and not in momentum space. The 
simplest model of spatial condensation that we know of is the model with a "defect" of 



Ref. [P^j . In this model M particles hop in a totally asymmetric way with a rate 1 on 
a ring with L sites and one particle (the "defect") hops with a rate p (with p < 1). It 
was shown that the partition function for the grand canonical ensemble is 



Z= E S(£n k -L)— exp U£ n * (CIO) 

n ,ni,... k=0 P \ k=0 ) 

where uq, ni, . . . are the number of vacancies in front of the defect, of particle 1, 
particle 2 etc. Obviously, 

L- (M + l) = ^-\ogZ (C.ll) 

At low density of particles (large density of vacancies), the particles produce a jam 
behind the defect and the string of vacancies (of length n ) in front of the defect can 
be seen as a condensate. The expression ( |C10| ) looks similar to ( |C.1| ) with at least two 
major differences: the occupations numbers n& are constrained by the length L of the 
chain (we have condensation in real space and not in momentum space!) and there are 
no measurable correlation functions given by creation and annihilation operators in this 
problem. 

In the model discussed in present paper, the analogy between the expression of 
the partition functions ( p. 32 ) for the stationary state and ( |C.1| ) for the Bose gas is 
completely obscure. Some similarities do exist. The algebra (p.2|) plays a role similar to 
the Heisenberg algebras hidden in the occupation numbers appearing in the Bose gas 
and the transformation (|2.4|) representing the conservation of the numbers of positive 
and negative particles (we have a U(l) x £7(1) symmetry in our case) corresponds to the 
transformation ( |C.3|) . We were not able to bring more understanding in the symmetry 
breaking mechanism using the algebraic approach although we wish we had, since we 
were looking for an equivalent of the Bogoliubov transformation in oder to eliminate 
the condensate. Our approach using Eq.( |9.6|) can be seen at this stage as an educated 
guess to do the regularisation. 
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